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Abstract 

Motivated by forthcoming p-Pb experiments at Large Hadron Collider which 
require both knowledge of gluon densities accounting for saturation and for processes 
at a wide range of pt we study basic momentum space evolution equations of high 
energy QCD factorization. Solutions of those equations might be used to form a set 
of gluon densities to calculate observables in generalized high energy factorization. 
Moreover in order to provide a framework for predictions for exclusive final states in 
p-Pb scattering with high pt we rewrite the equation for the high energy factorizable 
gluon density in a resummed form, similarly to what has been done in [T] for the BK 
equation. The resummed equation is then extended to account for colour coherence. 
This introduces an external scale to the evolution of the gluon density, and therefore 
makes it applicable in studies of final states. 



1 Introduction 

One of the scientific plans for the LHC is to collide protons with lead nuclei. This 
will offer an unprecedented possibility to test the partonic structure at high densities 
at relatively "cleaner" environment as compared to A A collisions [2]. The particularly 
interesting phenomenon which is believed to be understood with studies of observables 
produced in p-Pb is a phenomenon of gluon saturation [3]. For recent overview we refer 
the reader to [HE]. Already some time ago it has been recognized that saturation can 



be conveniently studied in the production of system of forward central di-jets [6HH], and 
recently that the signal of saturation might be expected in processes with high p t ~ 
20 GeV or larger [9]. With a system complicated enough like p-Pb, one may ask the 
question whether appropriate factorization of the process into hard parts, and parton 
densities can be establish taking all relevant degrees of freedom into account. To a good 
approximation i.e. neglecting the final state interactions of the forward jet one can use 
so called hybrid high energy factorization where one of parton densities parametrizes 
partons at large x while the other one at low x [TUHH]. However, it has been observed 
that generalized high energy factorization and appropriate hard matrix elements, and six 
gluon densities can be obtained providing more general framework. The six gluon densities 
as has been shown in [12] can be expressed in terms of two basic gluon densities which can 
be obtained from the BK equation [T3HT5] originally formulated in the coordinate space. 
Particularly useful for including higher order corrections, and numerical simulations of 
the scattering process is the momentum space representation/formulation of evolution 
equations. In this representation one can for example extend the BK equation to take 
into account sub leading pieces of the splitting function or bringing higher corrections 
in lnl/x kinematical effects [TBl[T7] . With momentum space formulation, one can, after 
performing resummation of the unresolved and virtual contributions in the kernel, extend 
the BK to include coherence which is a step towards studies of exclusive observables [1] . 

The resummation performed in [1] transforms the full kernel (linear, and nonlinear 
part) of the equation to a form where so called Regge form factor multiplies full ker- 
nel. This result is highly non-trivial since, the equation was nonlinear, and the resum- 
mation technique was based on an integral transform method usually applicable to lin- 
ear equations. The object for which the resummed equation was written was the so 
called Weizsacker- Williams gluon density (called in [1] dipole amplitude in the momen- 
tum space). This is one of the two basic densities which are used to construct densities 
used in studies of di-jets as derived in [T2"j . 

In the present article we investigate an equation of which the solution gives the so called 
high energy factorizable gluon density [18]. It turns out that although the nonlinear part 
of the equation for the unintegrated gluon density is much more complicated than the 
one studied in [1], one can perform still resummation analogously to [I], the results is the 
equation f l23|) . The new form of the equation under consideration allows for the extraction 
of a more compact expression for the triple pomeron vertex without splitting the vertex 
into real and virtual terms which will most probably be easier to handle in future numerical 
simulations. Although it is written in an exclusive form, obtained equation still does not 
allow for investigations of the exclusive processes since it does not include information on 
the hard process. The idea we follow is based on results from [T] where the BK equation 
was extended to include colour coherence by replacement of the Regge form factor with 
the non-Sudakov form factor, and the introduction of angular ordering. The new result for 
the unintegrated gluon density with angular ordering is given by equation f l28l) . For other 
approach to extend the equation for the unintegrated gluon density to be able to deal 
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with exclusive processes at high p t we refer the reader to [16] and for latest applications 
of this framework to [9]. 

The paper is organized as follows. In section 2 we recall the evolution equations for the 
high energy factorizable unintegrated gluon density, and the BK equation, and we show 
explicitly how to transform one to the other. In section 3 we perform resummation of 
virtual and unresolved real contributions in the linear part of equation for the high energy 
factorizable gluon density. Furthermore, we extend the obtained equation to account for 
colour coherence. We conclude in section 4. 

2 Nonlinear equation for high energy factorizable un- 
integrated gluon density 

In this section we will review and collect some facts about the equation for high energy 
factorizable unintegrated gluon density. We also provide some new insight on implications 
of its solution. The solution of the nonlinear evolution equation for the high energy 
factorizable unintegrated gluon density enters for example the formula for the proton 
structure function, which at leading a s \n(l/x) reads: 

F 2 (x,Q 2 ) = ^a s J2e 2 q J d 2 kF(x,k 2 ) (S L (k 2 ,Q 2 ,m 2 q ) + S T (k 2 ,Q 2 ,m 2 q )) (1) 
q 

The impact factors St, and Sl are obtained as matrix elements from the process 7*+g* — > 
qq with inclusion of appropriate phase space factors, the sum runs over quark flavours. 
In the massless case the full impact factor reads (see also [19] for massive case): 

s^ 0)+ s^ 0) - jf * j[ V - ^ - t^f^O 1K(1 = ° 

(2) 

For a full list of processes in which the gluon density J 7 is applicable we refer the 
reader to [T2] . 

In the formulas above Q 2 is the four momentum square of the photon, k 2 is the virtuality 
of the gluon, which satisfies k = |k|, I = |1| (where it will not lead to confusion we will 
not use for 2-d vectors bold font). Those momenta are lying in a plane transverse to the 
collision axis, (see Fig. ([1]) left). At the leading order in a s lnl/x the density T{x, k 2 ) 
obeys the following evolution equation [T6l[20l[2T] : 



/ dz rdi 2 



J r (x, k ) = Fq{x, k ) + a 
2a 2 n r 1 dz 



7 I I 2 

x/xq " JO h 



N C R 2 J x / X0 z 



"^(x/zj 2 ) 
Uk 2 1 



l 2 T(xlz,l 2 )-k 2 T(xjz, k 2 ) k 2 7{xjz, k 2 ) 
\k 2 -l 2 \ v /(4/ 4 + fc 4 ) 

r°° di 2 ( i 2 , 



F(x/z,k 2 ) I —\n\-)T(x/zJ 2 ) 



k 2 



(3) 
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Figure 1: Left: Kinematical variables. Right: Triple pomeron vertex convoluted with gluon 
density. 



where a s = N c a s /ir z — xjx' (see Fig. ([T]) left). The nonlinear term in ([3]) is a convolution 
of the angle-averaged triple pomeron vertex [20] with a gluon density. It reads: 



V{k, -k;h, -h,l 2 , ~h) 



™ 2 S 
N C R 2 



2Q(l\ - k 2 )B(l\ - k 2 ) + k 2 In - k 2 )6(l 2 2 - l{] 

+k 2 \n l h{il-k 2 )e{i 2 l -i 2 2 ) 

h 



(4) 



where the first term corresponds to the real contribution while the second and third to 
mixed virtual, and real terms. The arguments of the theta-functions indicate that the 
triple pomeron vertex orders the momenta in an anticollinear configuration. Convolving 
it with the gluon density (see Fig. (OQ) right), propagators l/l 2 , and l/l 2 (where remain- 
ing factors and l/l 2 were absorbed in the definition of gluon density), and changing 
variables we obtain the nonlinear term in equation ([3]) (for more detailed derivation of 
([3]) we refer the reader to [20] ). The strength of the nonlinear term is controlled by the 
parameter 1/R 2 which comes from the cut on a loop momentum circulating in the triple 
pomeron vertex when forward limit is taken. The inverse of a cut-off is interpreted as a 
hadron's radius [20J. This equation can be transformed to coordinate space by rewriting 
the convolution formula for the F 2 structure function as a convolution of the dipole am- 
plitude N(r,b,x) and the photon wave function [16]. Namely the relation between the 
dipole amplitude and high energy factorizable gluon density is [30] : 

F(x,k 2 ) = N " - [ d 2 b I d 2 re ik - r V 2 r N{r,b,x) (5) 



For later applications we also write explicitly the relation between the Waizecker- Williams 
gluon density and the dipole amplitude. It reads 

1 f d %f'ir 

2ir. I r 2 



$( X ,A; 2 ) = — I d 2 b I —e ik - r N(r,b,x) (6) 
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From the relations above it follows that 

J(x^ 2 ) = i 2 V^(x,fc 2 ) (7) 

which in principle could be defined up to terms of type: 

$ c = C 1 lnA; 2 //i 2 + C 2 (8) 

but the relation of k 2 ) to the dipole amplitude determines such terms to be zero. 
For more detailed discussion on this we relegate the Reader to the appendix. 
We also note that relation (J7J) is in fact a Poisson equation allowing to some extend for 
an interpreting analogously to electrodynamics J 7 as density while $ as a potential. The 
inverse transformation can be performed since (jSJ) does not contribute and the inverse is 
given by: 

a ir 2 r°° dl 2 I 2 

The function k 2 ) is linked to a function $&(a;, k 2 ) which obeys the BK equation (in 
momentum space): 



$ b {x,k 2 ) = <5>ob(x,k 2 ) + a s 



1 dz f°° dl 2 



x/x Z JO ^ 



l 2 <t> b {x/z, I 2 ) - k 2 $ b {x/z, k 2 ) k 2 $ b {x/z, k 2 )' 
\k 2 -l 2 \ + v /(4Z 4 + fc 4 T. 

-a s f —^ 2 b (x/z,k 2 ] 

(10) 



In the equation above, impact parameter is not a dynamical quantity. Furthermore (1101) 
does not describe physics at large distances correctly since it does not take into account 
confinement effects [23] . Taking above into account in order to obtain the integrated over 
h gluon density one has to assume some impact parameter dependence since one needs to 
perform the integration over it to calculate gluon density. For example assuming cylinder 
like target one may use the following anzatz for factorization [16] 



$ b (x,k 2 ) = ${x,k 2 )S(b). (11) 

with normalization conditions 



J d 2 bS(b) = l, J d 2 bS 2 (b) = -L (12) 



After integration over the impact parameter b we obtain 



$(x,k 2 ) = $ (x,k 2 ) + a s [ —f ^- 

Jx/XQ Z JO ' 



l 2 $(x/z, I 2 ) - k 2 <$>(x/z, k 2 ) k 2 <S>(x/z, k 2 ) 
\k 2 -l 2 \ + v /(4/ 4 + fc 4 J. 

dZ $ 2 (x/ Zj k 2 ). 

(13) 
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Figure 2: High energy factorizable gluon density. The maxima signalize emergence of the 
saturation scale Q s (x). Plot is made for x = 10~ 5 (continuous line), x = 10~ 6 (dashed line), 
x = 10 -7 (dotted line) 



where upon the integration over the impact parameter the inverse radius of the proton R 
multiplies the nonlinear term in (1131) . 

The transformation between ([3]) and ffl3|) is known, however it can be performed in a 
more straight forward way than in [23] and is presented below. The linear part of the 
kernel can be written in the form of a differential operator [15]. We have: 



k 2 ) = §o(x, k 2 )+a s 



dz 



-X 



x/x 



d 



3, 



log k 2 



$(x/z, k 2 )- 



a. 



irR 2 



dz 



<& 2 {x/z,k 2 ) (14) 



x/x 



where xil) = 2'0(1) — V'(t) — ^(1 — 7) is the characteristic function of the BFKL. We 
represent the BFKL kernel of the BK equation as a power series around 7 C = 0.373: 



X( 



d lnk 2)$(x,k 2 ) 



Xilc) + (-dink* ~ lc)x'(lc) + ^(-<9lnfc2 - 7c) Y'(7c) + - 



$(x, k 



(15) 

Inserting this expansion into f Tl^|) . and noticing that the BFKL kernel is invariant with 
respect to k 2 dk2k 2 d^ = df nk2 we obtain the linear part of the equation for the high energy 
factorizable unintegrated gluon density in the same form as for the WW gluon density. In 
order to invert the nonlinear part we use ([5]). Folding the linear kernel again we obtain: 



Fix, k 2 ) = Fo(x, k 2 ) + a s 



dz 

x/x z JO 



dl 2 \l 2 F{x/z, I 2 ) - k 2 F{x /z, k 2 ) k 2 T{x/z, k 2 ) 



I 2 



\k 2 



AN r R 2 k 



l 2 \ 

dz 



VWTk*) 



x/x 



k- 



n 2 



(16) 



where the nonlinear part is written in a new compact form. The equation above has quite 
interesting properties. In particular the main contribution to the nonlinear part comes 
from the anticollinear part of the triple pomeron vertex (recall the theta functions in the 
triple pomeron vertex). This has been already observed in [20J but the consequences of this 
finding were not fully clear. Since the nonlinear part enters the equation with a minus 
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sign, the anticollinear pole is subtracted from the linear kernel during evolution. The 
subtraction is stronger when the x is smaller where the importance of the nonlinear part 
of the equation is larger. The presence of the collinear and anticollinear contributions 
to the linear kernel of ( 1T6|) give rise to unordered emissions. The subtractions of the 
anticollinear pole via triple pomeron vertex result in the collinear ordering of the gluonic 
chain both in x, and in transversal momentum's most probable value which is given 
by emergent saturation scale Q s . This can be seen by inspecting the shape of the gluon 
density as function of k. The plot (J2J) shows that the lower values of x the saturation scale 
is harder and the distribution shifts towards higher values of the transversal momentum 
of a gluon with its maximal value at the saturation scale Q s . The property of providing 
most probable value is a feature that distinguishes the unintegrated gluon density obtained 
from (1T5|) from unintegrated distributions obtained from linear evolution equations or from 
(1T31) . The most probable momentum in the latter case is given by its boundary value, 
since the distribution diverges at small k, and some cut has to be applied. 

3 Resummed form of the equation for high energy 
factorizable gluon density 

In order to find the resummed form of ffT6j) we follow similar steps as in pQ. In principle 
as described above, (|T3|) and ( TT6l) are related via the Laplace operator so in principle after 
solving ( TT5j) one obtains (fT5|) . However, for practical future Monte Carlo simulation of 
p + Pb — > di—jets + X where both of these densities enter it is more convenient to have 
separate equations. Another reason for having both equations is that, the dynamics of 
gluon evolution is more transparent in ( JT6|) since one explicitly sees the role of the triple 
pomeron vertex in bringing the unitarity corrections via removing the anticollinear pole. 
Also the limiting cases like the double logarithmic limit is easily seen in the later. In the 
first step the angle dependence in the linear part ffT6j) has to be reintroduced, and virtual 
corrections have to be collected in a one expression. After this procedure one obtains an 
exclusive representation of the linear part, while the nonlinear we leave unchanged. Thus 
we have 

F(x,k 2 ) = F (x,k 2 ) + a s f 1 — J [^[T( x / z ,\k + q \2)-6(k 2 -q 2 )T(x/z,k)] 

Jx/xo Z \J 

•oo A1 2 p 

In —T[xjz,l 2 



AN r .R 2 k 



P k 2 ' 



(17) 



In the next step the linear part has to be decomposed into so called resolved, unresolved 
ones, and virtual emissions. The splitting is dictated by a singularity appearing in the 
denominator in the linear part of ( ITTl) when q — > 0. The parameter which controls this 
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splitting ji satisfies the following inequalities: q > ji for the resolved part , and q < \i for 
the unresolved part. Thus we have 



F(x, k 2 ) = J-o(x, k 2 

_ f 1 dz f d 2 q 



;i8) 



+ a s f - f ^(x/z, |k + q| 2 )% 2 - /i 2 ) 
Jx/x z J nq 

+ a s [ — [ ^| |k + q| 2 )#(/i 2 - q 2 ) - 6{k 2 - q 2 )T{xjz, k)] 

Jx/x z J nq 2L 



ira 2 k 2 
AN C R 2 



k 2 



-i 2 



I 2 



k 2 ' 



To perform resummation of virtual, and unresolved real emissions we are going to use the 
Mellin transform, defined as 



F(u, k 2 ) 

while the inverse transform reads 

F{x, k 2 ) = 



dxx J~ (Xj k 



c+ioo 



dcu x "^(oj, k 2 ) . 



Performing the Mellin transform, and using in the unresolved part |k + q| 
q 2 < ji 2 we obtain 



JF{u3, k 2 ) = jFq{oj, k z 
a s f (i 2 q 



(19) 

(20) 
k 2 since 

(21) 



+ ^ ^[^,|k + q| 2 )% 2 -^ 2 )] + ^ 



a s f d 2 q 



uj J q 



UJ 



J(co,k 2 )[e^ 2 -q 2 )-9(k 2 -q 2 )} 



f°° dV , V r , . 



in the nonlinear term we changed the variables x/z — >• y, and we integrated over x giving 
1/cj in front of the nonlinear part. We obtain 



JF(u), k 2 ) = J^oico, k l 

a s f d 2 (\— 



(22) 



+ — / |k + q| 2 )% 2 - /2 2 ) - -^(w, k 2 ) In - 



uj J 7rq^ 
4N C R 2 uj 



dyy w ~ l k 2 V 2 



uj 

00 dl 2 , I 2 



simplifying we obtain 



^(uj, k 2 ) = Pq{uj, k 2 



a. 



d 2 q 



uJ + uj4N c R 2 J 



uj + uj J nq z 
1 dyy^Vl 
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[^,|k + q| 2 )]% 2 -/i 2 



r di 2 , i 2 r . r2 



(23) 



where 



FnloJ, k 



ujFq{u, k 2 



bJ + U 

Inverting the transform we arrive at: 



k 2 

U = a s \n— . 

ft 2 



T(x, k 2 ) = JFq(x, k 2 



-a. 



dz 



A R (z,k,n)< 



x/xo 



—^—k 2 X7 2 



00 dl 2 

~~2 



Tiq 2 
I 2 



. x 



6(q 2 -fi 2 )^(-\k + q\ 

z 



(24) 



(25) 



-1 2 



where A R (z, k, jj) = exp ^— a s In - z In ^-J is the Regge form factor and we defined 



F (x, k 2 ) = — 



c+ioo 



dux~ w F°(uj,k 



From equation ( 1251) we may extract the new form of the triple pomeron vertex. 



rcsummcd 



2 12 

'2a 7„ ,.\V72t„ 'l Q(i2 ,2m„ l 2 /j/,2 ;„2^ 



Ri k 2 A R (z, k, In - k 2 ) In ^ - k 



(26) 



(27) 



In the form above it is understood that the Laplacian has to be taken after convolution 
with the parton density is performed. 

The last equation may be generalized to account for colour coherence, replacing the Regge 
form factor by the non-Sudakov form factor, and by introducing angular ordering in an 
angle of emitted gluons with respect to the beam direction. The extended equation takes 
the form: 



F(x, k 2 ,p) = F {x, k 2 ,p) 



(2? 



fd 2 q 

+ a 



f^f —0(p-qz)A ns (z,k,q){F(-,\k + q\ 2 ,q) 
J vq 2 J x/xo z [ z 



;q 2 S(q 2 



AN C R 2 

The momentum vector associated with i-th emitted gluon is 

q% = aip P + (3ip e + q ti 



(29) 



and thus we have rapidity and angle of emitted gluon with respect to incoming parent 
proton (beam direction) 



-In I — 

at; 



111 



\fsOLi 



8i 

tan — 
2 



\fso>i 



(30) 



The variable p in 028p is defined via £ = p 2 /(x 2 s) where |ln(£) is a maximum rapidity 
which is determined by the kinematics of the hard scattering, y/s is the total energy of 
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the collision. The particularly compact form of (128]) where the A ns acts both on the 
linear, and the nonlinear piece might be a subject to further extensions. For example 
one can use a non-Sudakov form factor which takes into account so called kinematical 
effects which bring important pieces of higher order corrections to those of ln(l/x) type. 
Another possibility is to account for sub leading pieces of the splitting function therefore 
proposing the extension of the CCFM equation [251428] for unintegrated gluon density 
with nonlinearities. 

The unintegrated gluon distributions obtained from equation ( 1231) corresponds to the 
coordinate space dipole in the fundamental representation. For LHC applications one also 
needs it in the adjoin representation. The corresponding unintegrated gluon density in 
the adjoin representation can be obtained by calculating the corresponding dipole cross 
section, and transforming it back to momentum space. To be precise one applies 

a(x, k 2 ) = ^j^(l- J (x, k 2 ))T(x, k 2 ) (31) 

together with 

N G (x, r, b) = 2N(x, r, b) - N 2 (x, r, b), a(x, r) = 2 J d 2 bN(x, r, b) (32) 

to obtain the gluon density in the colour adjoin representation: 

Fg{xX) = -£jnrs / d 2 b [ d 2 re^ k V 2 r N G (x,r,b) (33) 

After proposing equation (128]) one can generalize the formulas above to the situation 
where the dipole density depends on the hardness. In this case formula (13T]) changes to: 

8-7T 2 I" rlk 

a(x,r,p) = ^J ^(1- J (x,P))^(x,A; 2 ,p) (34) 

In a similar way equations below ( ]3T]) can be generalized to account for hardness. We 
have 

F G (x,k 2 ,p) = -%f- [ d 2 b [ d 2 re^ k V 2 r N G (x,r,p,b) (35) 
a s {2-K) 6 J J 

It will be interesting once the equation (]28j) is solved to check the behaviour of the dipole 
amplitude at large r, and in particular whether it flattens at constant value for different 
p. 



4 Conclusions, and outlook 

In this paper we studied evolution equations for the gluon density, the solutions of which 
form the basis for gluon densities entering the generalized high energy factorization frame- 
work. We have shown explicitly how the equations for the Weizsacker- Williams gluon 
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density, and the high energy factorizable gluon density are related to each other. In the 
next step we performed resummation of the unresolved real emissions in the evolution 
equation for high energy factorizable gluon density. The interesting aspect of our result is 
that the triple pomeron vertex itself does not need to be resummed since it is no singular, 
but the resummation of the linear part affects it in a simple multiplicative way. In the 
next step we proposed an equation which generalizes the resumed equation to account for 
colour coherence. The last equation is interesting because of possible future applications 
in studies of high p t exclusive final states with a focus on saturation. In the future it will 
be interesting to see via numerical calculations how the resummations performed affects 
both of investigated gluon densities i.e. how solutions of equations before resummation 
are related to solutions after resummation. On the level of linear equations it has been 
shown that the resummed equation gives for small enough fi identical answer to the solu- 
tion before resummation [21] and that the /x dependence is an effect of higher twist [28] . 
However, in the case of nonlinear equations the nonlinear part introduces effects of higher 
twists and therefore there might be some nontrivial interplay, for small enough saturation 
scale, of the saturation scale Q s and the cut-off scale fi. In will be also interesting to 
study to what extend the resummation affects the relations via Laplace operator between 
gluon densities i.e. what is the relation of direct solution of ( 12 5 p to action of laplacian on 
solution of resummed version of (TT3]) . 
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5 Appendix 

In this appendix we will perform the inversion of the Laplace operator in ([7]) and provide 
an example of action of it's inverse on the high energy factorizable gluon density within 
the GBW saturation model. We will also show how the behavior of the dipole amplitude 
sets constraint on the WW gluon density. The equation ([7]) reads: 

F( x ,k 2 ) = -^k 2 Vmx,k 2 ) (36) 

The physical boundary condition that the WW gluon density obeys in the highly per- 
turbative domain i.e. k 2 — > oo is that k 2 )\ k 2 =00 = and k 2 d^^(x, /c 2 )| fc 2 =00 = 0. It 
might be linked to the behaviour of the dipole cross section for small dipole size 

r r°° dr r r 1 ^ 2 dr 2 1 

$(x,k 2 ) = J d 2 b J -J (kr))N(x,r,b) = J d 2 b J —r 2 = S ± - (37) 
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where we used N(x,r,b) ~ 6 (bo — b) r 2 and J(0, kr) ~ 1 for small r [30l|35]. The theta 
function indicates that we model a hadron as a cylinder, the integral over the impact 
parameter b gives targets area S± while the r 2 behaviour of the dipole amplitude comes 
from colour transparency. The behaviour of the dipole amplitude at small r is sufficient to 
set the terms C\ and C2 to zero in (jSJ) since they would modify the gluon density globally 
i.e. the term In k 2 /fi 2 would lead to a divergent behaviour of the WW gluon density at 
large k while the C2 would introduce a constant contribution over the whole region which 
is also not consistent with 1/k 2 behaviour of the gluon density. Setting this terms to zero 
is therefore justified and allows for inversion of the Laplace operator [3lj with what we 
proceed below. The equation ( 136]) can be rewritten as 

F( X)U ) = -^- 2 d 2 u <$>{x )U ) (38) 

where u = \nk 2 jk\ we set k\ = 1 and denote A = • Integrating once and using the 
boundary condition <9 u $(x, u = 00) = we obtain: 

1 r°° 

d u $(x, u) — — - j du'F{x,u') (39) 

Integrating one more time and using the boundary condition for the function §(x,u = 
00) = we obtain: 

1 r°° 

$( X , M ) = _/ du'F{x,u'){u' -u) (40) 
Using the original variables we have 

Now we proceed with the example of the application of the above formula to high energy 
factorizable gluon density F(x, k 2 ). As an example we take the GBW gluon density [3T] . 
In this model the gluon density falls too fast at large k 2 as compared to perturabative 
QCD pattern but nevertheless it models quite well bulk properties of physics of saturation. 
It has also correct asymptotics and is frequently used as initial distributions for evolution 
equations P3JE3]. Applying the formula Q3SJ to F{x,k 2 ) = ^-k 2 /Q 2 {x)e- k ^ Q2 ^ we 
obtain 

Hx,k 2 ) = ^T(0,k 2 /Q 2 s ) (42) 

which satisfies the asymptotic behaviour $(x,/c 2 ) ~ for k 2 — > 00. The dipole cross 
section in the GBW model follows from inserting ( f42|) to inverse of ()6]) and reads 

r d 2 k 



a(x, r 2 ) = S ± r 2 / — e~ lk r ^(x, k ) = 2S ± ( 1 - ) (43) 



in agreement with [32] . 
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Alternatively to what we show above in order to determine the boundary conditions 
at large k 2 of the WW gluon density we can start with the lowest order expression in 
a s ln(l/rr) for high energy factorizable gluon density i.e. ^{x.k 2 ) = const-^. Using ( 13"8"|) 
we have 

~, ,9n a s^ 2 f°° dl 2 , I 2 ,o. a s n 2 1 , 1AS 

Hx - k) = -fcLT k k) = consi lvr^ < 44 > 

what is consistent with ( 1371) and ( 1421) and also justifies setting terms C\ and C2 to zero. 
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